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Abstract: We present results from numerical studies of the finite temperature phase tran- 
sition of the (3 + l)d 0(A^)-symmetric non-linear sigma model for = 1, 2 and 3. We study 
the dependence of the width of the 3(i critical region on A^ and we show that the broken 
phase scaling region is much wider for N = 2 and 3 than for A^ = 1. We also compare the 
widths of the critical region in the low T and high T phases of the 0(2) model and we show 
that the scaling region in the broken phase is much wider than in the symmetric phase. We 
also report results for the width of the scaling regions in the low T phase (2 + l)d Ising model 
and we show that the spatial correlation length has to be approximately twice the lattice 
temporal extent before the 2d scaling region is reached. 
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1. Introduction 



The behavior of symmetries at finite temperature is one of the most outstanding and relevant 
problems in current areas of particle physics such as cosmology, relativistic heavy-ion colli- 
sions and the quark-gluon plasma. It was shown some time ago that spontaneously broken 
symmetries in quantum field theories get restored at sufficiently high temperature In 
the high temperature limit, however, perturbation theory becomes unreliable when infrared 
divergences are not properly summed. Various groups developed temperature dependent 
renormalization group techniques to study the critical properties of the 0{N) scalar field 
theories (for a few examples see and references therein). 

In recent years considerable work has been done on the physics of the finite temperature 
chiral phase transition in QCD and related models. A compelling idea based on universality 
and dimensional reduction first put forward in ||^, is that in QCD with two massless quarks 
the physics near the transition can be described by the three-dimensional 0(4)-symmetric 
sigma model. Although, there is little disagreement that the chiral phase transition in QCD 
is second order, no quantitative work or simulations have been done that decisively determine 
its universality class, because simulations are performed away from the continuum limit on 
small volumes and show substantial discretization and finite size effects. The knowledge of 
the scaling relations are crucial for a successful extrapolation to small mass and large volume 
near T^. 

The reasoning behind the dimensional reduction scenario is based on the dominance of 
the light degrees of freedom near T^. In the imaginary time formalism of finite temperature 
field theory a given system is defined on a manifold x with the inverse temperature (3 = 

being the circumference of S"^. The boundary condition in the time direction is periodic 
for bosons and anti-periodic for fermions. At temperatures larger than the physical scales of 
the system, the nonzero Matsubara modes acquire a heavy mass and non-static configurations 
are strongly suppressed in the Boltzmann sum. If the original four- dimensional theory has 
fermions, then because of the anti-periodicity in the temporal direction, all fermionic modes 
become massive and can be integrated over. This is not to say that non-static modes 
have no effect. Rather they generate counter-terms to the three-dimensional theory. The 
physical picture of the dimensional reduction scenario can be described nicely by using 
a slightly different language presented in [Q. As one approaches the critical point, the 
fluctuations in space become correlated over a distance ^ which diverges at the critical point. 
In the critical region, not only do the fluctuations extend over large spatial scales, they also 
become very slow and there is a divergent relaxation time r; the effect is known as critical 
slowing down. The appropriate quantum of energy is then with ~ 1/t. Hence if the 
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transition is at a finite temperature Tc, hwc <S ksT near the critical point, which means 
that long wavelength fluctuations are classical thermal fluctuations and quantum mechanics 
is irrelevant. Equivalently, since the time interval (3 is very short compared to the relaxation 
time r, typical time histories will consist of a single configuration which is the same at 
each time slice. The dynamics therefore drops out of the problem and a Boltzmann weight 
exp(— /3ifciassicai) IS rccovcrcd from the path integral. 

In this work we study numerically the finite temperature transition both in the (3 + l)d 
and (2 + l)d 0^ model. This is equivalent to studying the phenomenon of dimensional 

crossover in a (i-dimensional layer 0{N) spin model of infinite extent in d — 1 dimensions 
and finite thickness in the remaining dimension. In other words temperature effects be- 
come predominant as a consequence of a change in dimension. The formal derivations of 
the dimensional reduction scenario require that the scalar mass is small compared to the 
temperature. We check whether this condition is necessary and we study the effects of the 
infrared fiuctuations on the width of the critical region under different circumstances. We 
performed simulations of the model with different 0{N) (for N — 1,2 and 3) symmetries 
in (3 + l)d and study the dependence of the width of the 3d scafing region on N. In the 
case of the (3 + l)d 0(2) model we compare the widths of the scaling region in the low 
temperature broken phase and the high temperature symmetric phase. We also study the 
scaling region of the (2 + l)d 0(1) = Z2 model. In addition we are interested to see how 
finite size effects can distort the results in such a simple model that can be simulated on 
lattices with relatively large temporal and spatial sizes. We extract critical exponents by 
measuring the order parameter, the susceptibility, and the correlation length. 



In order to study the problem of symmetry restoration in a non-perturbative way, we regard 
this model as an effective low energy theory by introducing a lattice regularization with a 
lattice constant a = 1/A, where A is the ultraviolet cut-off. Summing the bare Lagrangian 
over all lattice sites and transforming the continuum parameters into suitably defined lattice 
parameters according to a standard convention 



2. The model 



In the continuum limit the Lagrangian of the model reads: 




(2.1) 
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we are led to the following action: 



S|*| 



X 



(2.3) 



where the index /i runs over the nearest neighbors, (/x = 1, (i), and the index a runs over 
the components of the field $. For convenience, in our numerical work we concentrated 
on the limit of an infinitely strong bare coupling where the field $ is constrained to unit 
length: $ ■ $ = 1. In this case the action S reduces to the form: 



(2.4) 



X,fJ, 



which is equivalent to the energy of a configuration of the classical ferromagnetic 0{N) spin 
model. 



Finite physical temperature on the lattice means a finite size of the temporal extent Lt. 
The physical temperature is given by T = l/{Lta). The temperature of the system may be 
altered, by varying either the coupling k, which amounts to varying the lattice spacing or by 
varying Lt. To minimize the effects of lattice artifacts the condition A ^ T must be satisfied 
which requires a lattice temporal extent ^ 1. This condition not being met implies that 
the sum over the Matsubara modes is truncated at a small value. It is well-known that the 
{3 + l)d renormalized 0^ theory is a free theory. However, the regularized theory may well 
describe the correct physics of a broad spectrum of processes below the cut-off scale. 

As discussed in P] the regularized model can be interpreted in terms of physical terms by 
moving along the curves of constant physics (CCPs) which are characterized by g^i = const. 
Therefore, when approaching the T = second order critical line the mass of the sigma 
meson Mg- decreases and in order to keep T /M^ constant one has to increase Lt. Near the 
T = critical line the ratio T /Ma- should not change any more along a CCP. It was shown 
within a specific regularization scheme to leading order in 1/A^ [^] that Tc/M^ = \/SQ/{Ng{i) 
and by solving the gap equation on the lattice it was shown that the corrections due to 
a finite Lt are negligible for Lt > 4. In order to suppress possible discretization effects 
arising from 1/A^ corrections we performed our simulations with Lt = 12. The (2 + l)d (p^ 
theory has an interacting continuum limit and it has been shown that the non-linear sigma 
model is renormalizable in the 1/A^ expansion [0. The phenomenon of symmetry restoration 
at finite temperature in the three-dimensional theory was investigated by means of lattice 
simulations and finite size scaling techniques by Bimonte et al. . The authors of ^ located 
the intersections of critical lines corresponding to different values of Lt with various CCPs 
lying in the ordered phase of the T = theory and showed that for Lt > 3 the critical 
temperatures approach the finite physical temperature of the continuum theory. 
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Furthermore, the condition Ls ^ ^ must be satisfied in order to obtain results that 
approximate the thermodynamic hmit. For this reason at certain values of coupling k we 
performed simulations on lattices with different spatial extent in order to detect and control 
finite size effects. 

In our simulations we used the single- cluster algorithm because it reduces dramat- 
ically the correlations between successive configurations near a critical point. The cluster 
update is as follows: First choose randomly a refiection axis in the plane, denote the compo- 
nent of the spin $j that is parallel to this refiection axis by and that which is orthogonal 
by "l*^. Then choose randomly a site i of the lattice as a starting point for the cluster C. 
Visit all neighbor sites j of i and allow them to join the cluster with probability 

p{t,j) = 1 - exp[-K{<l>i<l>f + \<^i<^f\)]. (2.5) 

After this is done, visit all neighbors of the new sites in the cluster and add them to the 
cluster with probability p{i,j)- Iterate this step until no new sites enter the cluster. Finally 
fiip the signs of all $^ contained in the cluster. 



3. The observables 



In this section we present the various observables measured in the broken and symmetric 
phases of the model in our lattice simulations. The magnetization M" of a given configura- 
tion is defined as: 

M'^^yY.^:, (3.1) 

X 

where V = L^Lt is the space-time volume. However, in a finite volume and without the 
benefit of an explicit symmetry breaking interaction with an external magnetic field, the 
direction of symmetry breaking changes over the course of the run so that M" averages to 
zero over the ensemble. It is in this way that the absence of spontaneous symmetry breaking 
on a finite lattice is enforced. Another option is to introduce a field operator 



which is a projection of to the direction of M" separately for each configuration [|T0 
The effective order parameter 



^^l^{J2'^.,x) = {\M\), (3.3) 
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differs from the true order parameter extrapolated to the zero magnetic field limit by a 
factor ||lT|. One can easily show that the susceptibility of the order parameter can be 



expressed in terms of moments of the magnetization as follows: 

X = V{{\M\')-{\M\f). 



(3.4) 



In the symmetric phase one gets x = ^(I^P)- Cluster algorithms enable us to reduce the 
variance of the expectation values in the symmetric phase by using improved estimators 
[0, O. The improved estimator for the susceptibility is given by: 



Ximp 




(3.5) 



where \C\ denotes the number of spins in the cluster C. 

We determined the temporal exponential correlation length ^^^p of the sigma meson from 
the asymptotic exponential decay of the zero momentum connected two-point correlation 
function constructed from the field ^a,x, i-e., 



0(t) = 5^$., 



G{t) ^ (OoOi) - S2 A 



exp 



exp 



exp 



exp 



(3.6) 
(3.7) 



Similarly, we measured the spatial correlation length ^|^p from the exponential decay of 
G{x) in the spatial direction. The asymmetry of the lattice at nonzero temperature implies 
that ^Ixp 7^ ^exp5 because temperature effects on the mesonic self-energy become direction 
dependent. The breaking of Lorentz invariance implies that the spectral functions and the 
dispersion relations are modified in a non-trivial way by the thermal statistical ensemble. 
However, both ^^^^.p and ^g^^p are expected to diverge as T Tc according to the scahng 
law ~ (Tc — T)~'^ , where is a critical exponent of the corresponding 0{N) 3d spin model. 
One can construct correlation functions from various other operators 0{t) that couple to 
the scalar particle a, but according to [1^ the operators we defined in eq.( |3.6| ) give a better 
signal than other options. In the symmetric phase the two-point correlation function is 
G{t) = ((^i^{xi,t)) ■ ^j{^{xj,0)) \ and its improved estimator is given by 



G{t) 



imp 



N 

Jc\ 



5,,(C)$f<l>M, 



(3.8) 



where Sij{C) = 1 ii i and j belong to the same cluster otherwise Sij{C) = 0, [T^ 
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We also measured the so-called second moment correlation length ^2nd, defined by 



(3.9) 



where G{k) = 'Y2j{exp{ikxj)^Q^j) is the Fourier transform of the two-point correlation 
function. On the lattice ,^2nd is expressed as 



6 



f ix/F) 



1/2 



V4sin2(vr/L) 

where L = Lt or L^, F = G{k)\^^^^^^^^. The improved ^2nd 

requires (j(/c)imp 

given by 



(3.10) 



G{k) 




$^ sin(A;Xj 



(3.11) 



This estimator of the correlation length is very popular, because its numerical evaluation is 
simpler than that of the exponential correlation length which requires an identification of 
an asymptotic exponential decay. However, the two definitions of the correlation length do 
not coincide, since in ^exp only the ground state enters, while in the case of ^2nd according 
to eq.( |3.9| ) a mixture of the full spectrum is taken into account. As mentioned in |]TB[ the 
difference from one of the ratio ^exp/^2nd gives an idea of the density of the lowest states of 
the spectrum. It is well-known however, that in the high temperature phase the asymptotic 
exponential decay behavior of G{t) is reached at much smaller values of t than in the broken 
phase (see e.g. [0]). This behavior of G(t) implies that the difference between ^exp and ^2nd 
in the symmetric phase is smaller than in the broken phase. In our simulations we measured 
Gnd only in the symmetric phase and we observed that the difference between ^exp and ^2nd 
vanishes within error. 



4. Results 



In this section we present results from numerical simulations of the four-dimensional 0{N) 
{N = 1,2 and 3) non- linear sigma model at T 7^ 0. We also present results from simulations 
preformed in the low temperature phase of the three-dimensional Ising model. By measuring 
the correlation length ^ at the crossover into the 3d scaling region we study the dependence of 
the width of the critical region on the symmetry group 0{N) and on whether the symmetry 
is manifest (T > Tc) or broken (T < T^). We define the crossover coupling Across as the 
value of K that if it is included in the fitting window, the extracted exponent Praag does not 
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agree within statistical error with the exponent of the 3d model /3mag- we show in the 
next paragraphs, the change in the value of the exponent is accompanied by a substantial 
increase in the x^/d.o.f. of the fit, implying that the quality of the fit drops rapidly once we 
enter the crossover region. We managed to identify the value of Across with relatively good 
accuracy, because as we will see in the following paragraphs the value of k was varied by 
small increments at the crossover into the 3d scaling region. Throughout this work we used 
a large lattice temporal extent Lt = 12 in order to minimize discretization effects. 

First, we present the results for the (3 + l)d Ising model. For certain values of k, we 
monitored finite size effects on the order parameter by simulating the model on lattices 
with two different values of spatial extent Lg. For the final analysis we chose the data sets 
for which the systematic errors on E are less than 0.5%. More precisely for the range of 
couplings 0.15040 <k< 0.15100 we chose = 48, for 0.15010 <k< 0.15030 = 72 and 
for 0.14993 < n < 0.150075 = 96. For each value of n we typically generated 10^ sweeps 
of the algorithm. The magnetization data was fitted to the form 

S = kJ^-"^^. (4.1) 

The results for different ranges of n are presented in Table [I|. The simple power law fit is 
very good for data close to the transition, i.e. for k, < 0.15005. For these fits the exponent 
Pmag is consistent within error with a recent estimate of the critical index of the 3d Ising 
model /?j^^g = 0.3269(6) fl^. As shown in Table |l|, the value of /?mag deviates from the 
3d values and the quality of the fit gets worse as we include more data points at larger 
values of k. Therefore, we conclude that the crossover to the 3d universality class occurs at 
Across ~ 0.150075. The data and the fitting curve to the data with k < 0.150025 are shown in 
Fig. |1|. The statistical errors are smaller than the size of the symbols and for this reason we 
don't plot them in the figure. The data for the susceptibility and the correlation length for 
these values of the coupling were very noisy and did not allow us to measure the exponents 
7 and z/. 



range of n 






xVd.o.f. 


0.14993 - 0.150025 


0.32(3) 


0.149867(9) 


0.3 


0.14993 - 0.15005 


0.34(2) 


0.149863(7) 


0.5 


0.14993 - 0.150075 


0.37(1) 


0.149870(5) 


1.1 


0.14993 - 0.150100 


0.37(1) 


0.149853(5) 


7.3 



Table 1: Values of /3mag for the (3 + l)d Ising model extracted from different ranges of k. 



We measured the values of the spatial and temporal correlation lengths at the crossover 
into the 3d Ising scaling region by performing a simulation at k = 0.150075 on a 12 x 48^ 
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Figure 1: Order parameter vs. coupling for the (3 + l)d Ising model. The curve represents the fit 
to the data for k < 0.15005. 

lattice. For this measurement we chose a lattice with smaller spatial size Lg, because the 
correlation function was very noisy and we had to generate 15 million configurations in order 
to measure the correlation lengths ^* and C,'"^ with statistical errors less than 7%. By fitting the 
temporal and spatial correlation functions to the single-pole ansatz we extracted = 6.5(4) 
and = 8.8(5). The overall error was estimated from the observation of the fit sensitivity to 
the variation of the parameters and to the range of points used in the fit. We conclude that 
in order to enter the 3d scaling region it was not necessary to work so close to the transition 
itself that the correlation lengths would be large compared to the temporal extent of the 
lattice. The importance of this comment lies in the fact that the formal derivations of the 
dimensional reduction scenario require that the zero mode dominates the Matsubara sum 
and this occurs when the scalar mass is small compared to the temperature. Apparently, 
our result shows that the higher Matsubara frequencies do not affect the critical singularities 
but they do affect non-universal aspects of the transition such as the width of the scaling 
region. 

As the temperature approaches Tc the mass of the sigma meson decreases to zero and 
this results in an increase of the mcsonic number density. A consequence of this is that 
the spectral function may become broader due to collisional broadening. One can study 
in detail the properties of the mesons at T 7^ by using the maximum entropy method 
to extract their spectral functions from high statistics correlation functions measured on 
anisotropic lattices with temporal sizes larger than Lt = 12. However, for our purposes 
the Lf — 12 temporal extent is large enough for an estimate of ^* at Across- To the extent 
that we recognized particle-like excitations at T 7^ and ^ imphes that at Across the 
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breaking of Lorentz invariance caused a significant "renormalization" of the speed of light 
in the mesonic dispersion relation. 



range of k 


Ls 


/^mag 




xVd.o.f 


0.3125 - 0.3275 


24 


0.340(1) 


0.30358(7) 


0.5 


0.3125 - 0.3275 


48 


0.343(2) 


0.30362(1) 


1.1 


0.3125 - 0.3325 


48 


0.335(1) 


0.30390(6) 


5.6 


0.3125 - 0.3375 


48 


0.329(1) 


0.30417(4) 


7.8 



Table 2: Values of /9mag for the (3 + l)d 0(2) model extracted from different ranges of n. 
0.6 




Figure 2: Order parameter vs. coupling for the (3 + \)d 0(2) model measured on 12 X 24^ and 
12 X 48^ lattices. The curve represents the power law fit to the data for k < 0.3275. 

In this paragraph we discuss the numerical results in the low temperature phase of the 
0(2) symmetric sigma model. We generated approximately 80,000 configurations for each 
value of the coupling on the 12 x 24^^ lattice and approximately 40, 000 configurations for 
each value of k on the 12 x 48^ lattice. The results for the fits of S to the standard scaling 
law for different values of k are shown in Table |^. By monitoring the change in the quality 
of the fits and the deviation in the value of /^mag from recent estimates of /5^^g = 0.3485(2) 
[|15|] we conclude that Across ~ 0.3325. An interesting observation is that finite size effects 
are negligible even when the asymmetry ratio is two. In Fig. |^ we plot the order parameter 
versus coupling. The continuous line represents the power law fit to the 12 x 48'^ data for 
K < 0.3275 . 



The values of the spatial and temporal correlation lengths measured at the crossover 
coupling Across = 0.3325 are Across = 0.94(2) and Across — 0.97(3). We conclude that in the 
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0(2) model the crossover to the 3d universahty class occurs at a much smaller value of 
C,/Lt than in the 0(1) model. This explains why the critical region in this case is more 
easily accessible on lattices with a smaller asymmetry ratio than in the case of the Ising 
model. The fact that the value of Across is very close to the value of Across implies that 
the thermal fluctuations are relatively small at the crossover into the 3d scaling region. 
The diverging correlation length of the Goldstone boson is a source of additional finite size 
effects in a system. In this case however, the values of the order parameter in the 3d scaling 
region measured on small lattices are close to its thermodynamic limit values, because the 
Goldstone boson had the effect of increasing the width of the critical region. 
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Figure 3: Susceptibility vs. coupling for the (3 + l)d 0{2) model measured on 12 X 24^ 12 X 48^ 
and 12 x 72^ lattices. 



In Fig. ^ we present the results for the susceptibility x of the order parameter for the 
0(2) model. As expected finite size effects and statistical errors are much bigger for x than 
for S and for this reason we are unable to extract a meaningful value for the exponent 7. 



range of n 


Ls 


/^mag 




xVd.o.f 


0.4725 - 0.4875 


48 


0.377(4) 


0.4569(2) 


2.5 


0.4725 - 0.4875 


72 


0.364(4) 


0.4574(2) 


0.4 


0.4725 - 0.4975 


72 


0.356(2) 


0.4579(1) 


3.5 


0.4725 - 0.4500 


72 


0.351(1) 


0.4580(1) 


4.0 



Table 3: Values of (3mag for the (3 + l)d 0(3) model extracted from different ranges of k. 



We also performed simulations of the 0(3) model on lattices with sizes 12 x 48^ and 
12 X 72^ and we observed a similar scenario to the one described above for the 0(2) model. 
Again, we generated approximately 30, 000 to 50, 000 configurations for each value of k 



-10- 




0.45 0.46 0.47 0.48 0.49 0.5 0.51 0.52 0.53 0.54 0.55 



Figure 4: Order parameter vs. coupling for the (3 + \)d 0(3) model measured on 12 X 48^ and 
12 X 72^ lattices. The curve represents the power law fit to the data for k < 0.4875. 

on both lattice sizes. The values of /9mag extracted from fits near the transition are very 
close to recent numerical results which predicted /3^g = 0.3689(3) ||16[. The discrepancy in 
the values of /3mag extracted from the simulations on the two different lattice sizes and the 
difference in the value of x^/d.o.f. could be attributed to finite size effects on the smaller 
lattice. As expected the finite size effects in the 0(3) model are larger than in the 0(2) 
model, because it has twice as many Goldstone bosons. Again by monitoring the quality of 
the fit and the deviation of Pmag from /5^g we found that the crossover into the 3d scaling 
region occurs at Across ~ 0.4975. The values of the spatial and temporal correlation lengths 
measured on a 12 x 24^ lattice at the crossover coupling Across = 0.4975 are Across = 0.92(3) 
and Across = 0.90(2). Given the uncertainty in the measurement of the Across, we conclude 
that the crossover in the 0(2) and 0(3) models occurs at approximately the same value of 
the correlation length. 



range of n 


Ls 


statistics 


0.22425 - 0.22800 


240 


100,000 


0.22404 - 0.22417 


480 


300,000 


0.22400, 0.22404 


720 


500,000 



Table 4: Lattice sizes and statistics at different values of n for the (2 + l)d Ising model simulations. 

We also performed simulations of the (2 + l)d Ising model. The results for the order 
parameter are shown in Fig. ^. In this case, because the dimensionality is lower, the infrared 
fluctuations are stronger than in the (3 + l)d model with the same symmetry. Therefore, in 
order to suppress finite size effects we simulated the model on lattices with larger asymmetry 
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ratios than in the (3 + l)d case. The values of Lg and tlie amount of statistics generated at 
various values of n are shown in Table ^ By fitting the data with the four smallest values of 
K to the standard scaling relation we extracted Kc = 0.22394(1) and /3mag = 0.138(8), which 
is close to the 2d Ising exponent /5^^g = 0.125. We conclude that the crossover into the 
2d scaling region is at Across ~ 0.22417. The value of the spatial correlation length at Across 
measured on a 12 x 240^ lattice by generating 5.5 million configurations is Across = 25.0(2.0). 
Unfortunately, the large statistical fluctuations in the data set did not allow us to measure 
the temporal correlation length. 

Next, we present results in the symmetric phase of the 0(2) model. We measured the 
improved estimators for the susceptibility (see eq. ( p.5|) ), for the exponential correlation 
length ^exp extracted from the asymptotic decay of C(t)iinp (see eq. (p.8|)) and for the second 
moment correlation length ^2nd (see eq. ( ^.1U| )). It is well-known that in the high temperature 
phase of this model Ls/C, ~ 7 is sufficient to give thermodynamic limit results. In order to 
detect finite size effects, we performed simulations at /t = 0.2930 on lattices with = 24 
and 36 and we found that in both cases the values of the susceptibility and the correlation 
lengths agree within errors, implying that finite-size effects are under control. In most cases 
we kept Ls/^ > 9 by increasing Lg from 24 to 96 as we approached the critical coupling 
with constant Lt = 12. For each value of the coupling we generated approximately 10^ 
configurations. We also checked that the deviations between ^exp and ^2nd are less than 0.5% 
for all the values of k. An interesting observation is that ^* and remain equal for all values 
of K studied in this work. The value of k which is closest to Kc is 0.3015 and ^ = 14.18(7) 
at this value of the coupling. This observation implies that in the symmetric phase the 
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quantum fluctuations remain more important than the thermal fluctuations up to relatively 
large values of ^/Lt. 



range of k, 


range of 


7 


V 


0.2800 - 0.2970 


1.691(2) - 4.003(7) 


1.11(1) 


0.602(4) 


0.2990 - 0.3015 


5.33(3) - 14.18(7) 


1.18(1) 


0.59(1) 



Table 5: Values of 7 and v measured in the symmetric phase of the (3+ l)(i 0(2) model in different 
ranges of k and 




0.28 0.285 0.29 0.295 0.3 0.305 

K 



Figure 6: Correlation length vs. coupling measured in the symmetric phase of the (3 + \)d 0(2) 
model. The continuous line represents the fit to the data for 0.2800 < k < 0.2970 and the dashed 
line represents the fit for 0.2990 <k< 0.3015. 

We also measured the exponents 7 and v in the 0(2)-symmetric phase by fitting our 
data to the standard scaling laws 

X = B{k, - kY^- i = C{k, - (4.2) 

The results for the fits in different scaling windows are shown in Table |^. The data and the 
fitting curves are shown in Figs. ^ and ^ The exponents measured are significantly different 
from recent predictions of the three-dimensional critical exponents 7^^^ = 1.3177(5) and 
v^"^ = 0.67155(27) Our results are consistent however, with the observation mentioned 
above that in the symmetric phase the thermal fluctuations do not influence decisively the 
critical properties of the model even up to values of the correlation length that are much 
larger than the value of Across in the broken phase. 
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Figure 7: Susceptibility vs. coupling measured in the symmetric phase of the (3 + l)d 0{2) 
model. The continuous line represents the fit to the data for 0.2800 < k < 0.2970 and the dashed 
line represents the fit for 0.2990 < n < 0.3015. 



5. Conclusions 



We presented results from numerical simulations of the 0{N) sigma model at finite temper- 
ature for A?" = 1, 2, 3 in four dimensions and N — 1 in three dimensions. We investigated 
the phenomenon of dimensional reduction and studied the dependence of the width of the 
critical region on the number of field components N and on whether the 0{N) symmetry is 
spontaneously broken or restored. Contrary to the naive geometric scenario which requires 
that the condition ^ ^ Lt must be satisfied in order to reach the d — 1 scaling region, we 
showed that in the low temperature phase the correlation length of the sigma meson doesn't 
have to be much larger than Lt in order to reach the dimensionally reduced critical region. 
Our results imply that the nonzero Matsubara modes do not affect the critical singularities 
but they can affect non-universal aspects of the transition such as the width of the scaling 
region. We measured the correlation length of the sigma meson in the temporal and spatial 
directions at the crossover into the d—1 scahng region and we found that lor ^ — Lt — V2 
in the case of the 0(2) model {*j,oss = 0.94(2) and {^^oss = 0.97(3) and for the 0(3) model 
'Across ~ 0.90(2) and ^^ross ~ 0.92(3). We should note however, that the deviations from the 3(i 
scaling laws at these relatively small correlation lengths should also be affected by discretiza- 
tion effects. In the case of the Ising model wc measured ^^ogg = 6.5(4) and .^^ross ~ 8.8(5), 
implying that the critical region is much wider in the presence of the infrared fluctuations 
due to the Goldstone bosons. The critical regions for = 2 and 3 have approximately the 
same width and therefore we don't expect a significant change in the width of the scaling 
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region for > 3. The fact that for the continuous symmetries Across ~ Across imphes that the 
thermal fluctuations are still very small at the crossover temperature. 

It is not yet known whether the dimensional reduction scenario is the correct description 
of the QCD thermal phase transition. If it is, then the thermal transition with massless 
quarks may belong to the 0(4) universality class, whereas when the quarks are massive 
the tricritical point in the (/i, T) plane becomes critical and it must belong to the 3d Ising 
universality class. According to our results the latter must have a much thinner scaling 
region than the massless quarks transition. 

We also showed that the critical region of the (3 + l)d 0(2) model in the symmetric 
phase is much thinner than in the broken phase. Despite the fact that we simulated the 
model very close to the critical point and went up to correlation lengths ~ 14, the values of 
the measured exponents are not consistent with the 3d exponents; they are rather effective 
exponents with values between mean field exponents and 3d exponents. 

We also studied the symmetry restoration transition in the (2 + l)d Ising model. This 
model has a well-defined interacting continuum limit and its behavior near is particularly 
hard to explore analytically, because infrared divergences get worsened in three dimensions. 
We showed that in order to reach the 2d scaling region on lattices with temporal extent 
Lt = 12 the spatial correlation length C,^ = 25.0(2.0) and in order to minimize finite size 
effects the simulations near were performed on lattices with large spatial extents that 
are larger by almost an order of magnitude than in the four-dimensional case. In this lower 
dimensionality we studied only = 1, because according to the Coleman- Mermin- Wagner 
theorem for any continuous symmetry in two spatial dimensions the infrared divergent 



fluctuations wash out the order parameter and the symmetry is manifest for T > 0. 

In the future we would like to study the width of the scaling region of the 0{N) linear 
sigma model and compare with the results from the nonlinear sigma model presented in 
this paper. It was shown in the past that the dimensional reduction scenario is the correct 



description of the chiral phase transition in four-fermion models [0, because the infrared 
region is dominated by the zero Matsubara mode of the bosonic field and the nonzero bosonic 
and fermionic modes decouple. Simulations of the (3 + l)d four-fermion model with Z2 and 
SU{2) X SU{2) chiral symmetries are underway in order to study the dependence of the 
width of the critical region on the degree of chiral symmetry of the model. 
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